Yang-Lee zeros and the critical behavior of the infinite-range two- and three-state 

Potts models 
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The phase diagram of the two- and three-state Potts model with infinite-range interactions, in 
the external field is analyzed by studying the partition function zeros in the complex field plane. 
The tricritical point of the three-state model is observed as the approach of the zeros to the real axis 
at the nonzero field value. Different regimes, involving several first- and second-order transitions of 
the complicated phase diagram of the three state model are identified from the scaling properties 
of the zeros closest to the real axis. The critical exponents related to the tricritical point and the 
Yang-Lee edge singularity are well reproduced. Calculations are extended to the negative fields, 
where the exact implicit expression for the transition line is derived. 
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I. INTRODUCTION 

Since the famous work of Yang and Lee [l], Q , the zeros 
of the partition function in a complex activity plane were 
extensively studied as a useful tool in the investigation of 
various aspects of phase transitions, up to giving rise to 
a new type of criticality related to complex parameters 
such as the Yang-Lee edge singularity Q. Although orig- 
inally introduced in context of Ising model for explaining 
singularities related to the second-order phase transition, 
zeros appeare to be useful for deriving the related critical 
exponents and even distinguishing the nature of a phase 
transition for a variety of models (see e.g. [H, or Q for 
a recent review) . 

In the case of the Ising model, where all the zeros lie 
along the circle in the complex plane, the correspondence 
between the density of zeros and thermodynamic poten- 
tials is simple and straightforward. For the Potts model 
Q with an arbitrary number of states, the layout of Yang- 
Lee zeros is much more complicated |7Hllj while even 
more intensive investigations were conducted on Fisher 
zeros [l2l . [l3| in this model [l4| and the zeros in the 
complex-g plane [l5l - fl8l |. As far as the Yang-Lee zeros 
are concerned, it was shown analytically in one dimen- 
sion fl9| , and by numerical calculations on small systems 
in higher dimensions [|| [Icj, that, for general q, these 
zeros do not lie on the unit circle of the complex activ- 
ity plane, and that multiple lines, or bifurcations may 
occur This is related to a more complicated phase 
diagram of this model. When the number of states of 
the Potts model is larger than 2, the presence of real, 
positive external field h does not necessarily destroy the 
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phase transition, but produces a line of first-order phase 
transitions ending with the tricritical point. The model 
is not symmetric under exchange h -H- —h, and in the 
presence of a negative field, a different and less explored 
|2fj | transition takes place. All this leads to a more com- 
plicated phase diagram in the temperature-field plane, 
that involves both continuous and first-order phase tran- 
sitions, and the tricritical points. In the complex activity 
plane, this is reflected by the fact that zeros do not lie 
on the unit circle, while more than one line of zeros that 
approach the real axis may appear. 

The main purpose of this paper is to present, in the 
special case of the three-state Potts model with infinite- 
range interactions, an integrated view covering different 
critical regimes of the entire phase diagram in terms of 
the behavior of the partition function zeros in the com- 
plex activity plane. We also report the results for the 
Yang-Lee edge singularity above the tricritical tempera- 
ture of this model. 

The infinite-range Potts model is equivalent to the 
mean-field (MF) approximation of the standard Potts 
model with short-range (SR) interactions dHHH, which 
describes the critical beahaviour of the SR model cor- 
rectly only in sufficiently high dimensions, above the 
upper critical dimension. Nevertheless, its phase dia- 
gram displays similar diversity to that of the SR model, 
while it may be approached analytically, using the saddle- 
point approximation. In our approach we apply twice the 
Hubbard-Stratonovich transformation to obtain the ex- 
act two-parameter analytic expression for the free energy 
of the infinite-range model. The same approach permits 
us to obtain implicit analytic expressions for the parti- 
tion function zeros and to calculate zeros close to the real 
axis in the finite system. The finite-size scaling analysis 
of a few zeros that lie closest to the real axis then allows 
us to identify different critical regimes, and calculate the 
related critical exponents. 
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The outline of the paper is following. In Sec. |TT] we 
illustrate the approach for the simple q = 2 case, which 
is also used to obtain some estimates about the precision 
of numerical calculations of the partition function zeros. 
The same methodology is then applied to the q = 3 model 
in Sec. Mil We first derive an exact two-parameter ex- 
pression for the free-energy density to be solved by a 
saddle-point method, and we discuss its shapes in the 
two-parameter space, related to different regimes of the 
phase diagram. In Sees. MI Al and IIIIB1 we derive ana- 
lytical expressions for the scaling properties of the lowest 
zeros in different critical regimes, for the cases with h > 
and h < respectively. In Sec. MIB II we present the nu- 
merical results for the partition function zeros, and the 
numerical derivation of scaling exponents. Section IIVI 
contains the concluding remarks. 



TABLE I. Exact values of anomalous dimensions yu and y g 
as a function of K. 
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It may be solved by a standard saddle-point approxi- 
mation, which we recall here briefly. The corresponding 
extremal condition may be expressed by 



S = tanh 



KS + h 



(9) 



II. q = 2 (ISING) MODEL 

The energy of the g-state MF Potts model with TV 
particles in the external field H has the form 
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where Sj denotes the particle at site i, which can occupy 
one of the q Potts states, while any two particles interact 
with the same two-particle interaction Jo /TV, regardless 
of their distance. 

The reduced energy may be written in terms of num- 
bers of particles belonging to each of the q states, which 
for q = 2 are denoted by TVi and TV 2 
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and, by conservation of the total number of particles, 
TV = TVi + TV 2 . The partition function 
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where S = (2 TVi/TV) — 1 is the order parameter, and the 
model has a second-order phase transition for h = and 
A' = 2. 

The density of partition function zeros, g, which all are 
located on the imaginary axis of the complex field plane, 
is related (23Tj25| to the jump in the order parameter S, 



S = 2ng, 
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where S is the solution of Eq. (|9|). Above the critical 
temperature (K < 2) the gap free of zeros opens around 
the real axis. At its edge /i2,edoe; which depends on tem- 
perature as Q 
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another second-order phase transition in this model takes 
place, known as the Yang-Lee edge singularity. 

The exact expressions for the critical exponents are 
known 0, H2| for both transitions. We shall be inter- 
ested in scaling properties of the field h and the density 
of zeros g for which the exact values of their anomalous 
dimensions in different critical regimes are summarized 
in Table |TJ Note that the anomalous dimensions y~h and 
y g are not defined here in the standard way. Due to the 
infinite range of interactions the space scale has no mean- 
ing, so that the size scaling is carried out by scaling the 
number of particles instead. To to obtain the standard 
anomalous dimensions for this model as obtained when 
applying the mean-field approximation to the short-range 
version of this model, one should multiply fl8| the above 
exponents by d c , the upper critical dimension of the cor- 
responding phase transition (which is equal to 4 for the 
Ising transition and 6 for the Yang-Lee edge singularity). 

The anomalous dimension of the field is denoted by 
Uh, which may be obtained from the scaling with size of 
the distance (in the parameter space) of the imaginary 
component of the field hi from its edge value, h e( i ge 
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FIG. 1. (Color online) Solutions of the equations ReZN = 
[blue (dark gray) lines] and ItuZn = [orange (light gray) 
lines, including the hi = axis] for N — 1 000 000 particles, at 
K — 2. For distant points errors can be observed as deviations 
of intersections from the imaginary axis. 



The exponent y g = y^a denotes here the anomalous di- 
mension of the density of zeros, the power-law decay of 
which is described near criticality by the exponent a 
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Let us now sketch briefly the numerical derivation of 
these scaling properties by using the partition function 
zeros closest to the real axis. 



A. Partition function zeros - Numerical results 

For finite N, the partition function is the polyno- 
mial of order N in the complex variable z = e . The 
usual root-finding routines soon become inapplicable for 
finding the roots of Z^, as iV increases. Instead, we look 
for the numerical solutions of the two equations 
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in a complex (hi, hi) plane. A typical example of layout 
of these zeros, given by the intersections of dark gray 
lines and light gray lines hi = 0, is presented in Fig. [T] 

Consistent with the Yang-Lee theorem, these zeros lie 
on the imaginary axis. As can be observed in Fig. [TJ the 
precision of such a calculation deteriorates for the remote 
zeros, which start to deviate from the imaginary axis, Q, 
as a consequence of the accumulation of numerical errors, 
which increase with N and with decreasing K. However, 
we do not have to deal with this problem here, since 
only a few zeros closest to the real axis suffice to detect 
different regimes and to analyze the corresponding singu- 
larities. Our analysis will thus be focused on these zeros, 
and include a comparison to the exact results allowing 
the estimation of related errors. 

We calculated the position of the two zeros closest to 
the real axis for a sequence of different sizes, ranging from 



TABLE II. Numerical (/i2,extr), and analytical ( /i2,ed g e) values 
for the Yang- Lee edge, followed by the numerical convergence 
exponent j/2 and its exact value yu from Table [T] The error 
bars are estimated to be of order of the last cited digit. The 
numerical part of the errors is at least one order of magnitude 
smaller than the overall errors (including also the finite-size 
corrections) . 



K 


/l2,extr 


^2, edge 


2/2 


Vh 


2.5 


10" 9 





1.00002 


1 


2.0 


10~ 8 





0.7497 


3/4 


1.9 


0.01514 


0.015136 . . . 


0.664 


2/3 


1.5 


0.1814 


0.18117... 


0.65996 


2/3 


1.1 


0.478 


0.4756 . . . 


0.67 


2/3 



10 2 up to 10 6 particles. The distance of the zero closest 
to the real axis is used as an estimate of the position 
of the Yang-Lee edge at a given temperature, while its 
convergence to its value in the thermodynamic limit is 
expected, according to scaling property (fl"2j) . to give the 
anomalous dimension y^ 
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The results of extrapolating data by a simple least- 
squares fit are presented in Table ITTl 

The errors of the extrapolated values in Table [TT] (as 
well as in all the following tables) have two origins. First 
is the pure numerical error, which is relatively small (e.g., 
for 7/2 it ranges from 10~ 4 to 10~ 6 with increasing K) due 
to the high accuracy (10 — 12 digits) of the input data. 
The second type of errors is due to the finite-size effects, 
and it comes from neglecting the higher order correction 
terms in (| 1 5 1) and f) 18[) . They are in all cases one or more 
orders of magnitude larger than the numerical errors. In 
the tables, the extrapolated values are presented up to 
the first digit that differs from the exact value. Results 
presented in the second column clearly indicate the open- 
ing of the gap for K < Kq = 2. Good agreement also 
exists between the convergence exponent y 2 and the ex- 
actly known anomalous dimension y^ of the field (fourth 
and fifth columns). 

The second quantity of interest is the edge density of 
zeros, qn, defined as 
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The density of zeros near the Yang-Lee edge can be nu- 
merically calculated from the distance between the two 
zeros closest to the real axis, h^ and h^K 
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using a discrete version of Eq. (|16|) with dn = 1. It was 
extrapolated to the limit N — > oo by a least-squares fit 



TABLE III. Numerical (g C xtr) and analytical [(/exact, from Eq. 
(|10p ] values for the edge density of zeros, followed by the 
numerical exponent j/ 9 , e xtr and its exact value y g (from Table 
U). The error bars are estimated to be of order of the last cited 
digit. The numerical part of the errors is at least one order of 
magnitude smaller than the overall errors (including also the 
finite-size corrections). 
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to the form given by Eq. (|T3|) . i.e 
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The results of these extrapolations are presented in 
Table Mil The vanishingly small numerical values of the 
edge density of zeros indicate a second-order transition 
for K < 2. At lower temperatures, the transition is of 
first order, and the edge density of zeros has a finite value. 
The critical exponent of the Yang-Lee edge singularity is 
in excellent agreement with the exact value for the Ising 
phase transition, given by the scaling at K = 2. The 
scaling exponent for K < 2, corresponding to the Yang- 
Lee edge singularity, indicates a different critical regime, 
but is obtained with less precision. 



III. 



3 MODEL 



For the three-state Potts model in the symmetry break- 
ing field, the reduced energy may be written as 
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N 1 +N 2 + N 3 = N, 
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where Ni,N 2 , A3 denote numbers of particles in the three 
respective states. The external field is conjugated here 
to the Potts state 1 and (dis)favors it, depending on the 
sign of h. 

Compared to the Ising case, this model has a much 
richer phase diagram (Fig. [5]), the details of which will be 
discussed later in this section. In addition to the number 
of exact results for h > and arbitrary value of q, [22l |. 
we pay special attention to the less investigated h < 
part of the diagram. 

For h = and q > 2, the Potts model has a first-order 
phase transition at inverse temperature 
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FIG. 2. (Color online) Phase diagram: Ising transition 
(h < 0), open circles; KMS transition [h > 0), thick full 
line. At low temperatures, K > Ko(3), Ising and KMS lines 
are merged and a transition exists only for h — 0, broken line. 



For g = 3, K„(3) = 4 In 2 = 2.772588722239781 . . .. 

While in the presence of an external magnetic field, 
there is no transition for q = 2, models with q > 2 conceal 
a more complicated critical behavior, [26h29] . with two 
types of field-driven transitions: for h > and for h < 0. 
For h > 0, there is a line of first-order transitions in 
(K, h) plane 
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(22) 

which starts at (Ko(q), 0) and ends at the tricritical point 
(K tr , htr) with 



K tr (q) = 4- 
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-, (23) 

where the transition is of second order. For q = 3, 
K tr (3) = 8/3 = 2.6, h tr (3) = 0.02648051389327864... 
and the transition line is given by thick full line in Fig. 
® 

For negative values of h, the line of first-order transitions, 
also starts at (Ko(q),0) and approaches asymptotically 
to the point (Ko(q—l), —00) (open circles in Fig. [5]). The 
exact functional dependence h = h(K) is not known. 

For the purpose of numerical calculations, it is useful 
to present the partition function in a polynomial form 
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On the other hand, for the analytical approach, the par- 
tition function, written as 



e K(N-l)/2 e -(N K/4)[(4/3)(h/K-l/2) 2 + l] 
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allows a twofold application of the Hubbard-Stratonovich 
transformation ([5]) which reduces it to a double integral 
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with the exact free-energy density f(x, S) 
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where S = (3Ni/N — l)/2 is the order parameter for 
h > 0, and x = (N 2 — N^)/N is the order parameter for 
h < 0. The integral (|27j) is then solved by a saddle-point 
approach. The locations of minima are the solutions of a 
two-by-two system of equations for x and S 
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The MF approaches to the general q-state Potts model, 
such as the solution by Kihara et a/.[2l[ usually neglect 
the fluctuations among states orthogonal to the ordered 
one, which corresponds to taking x = 0. We examine 
here the free energy in the entire (x, S) plane. In Figs. [3] 
- 16] we illustrate the shapes of / for several characteris- 
tic values of temperature and the field, corresponding to 
different phases. For clarity in the figures, the maxima 
of — / are displayed instead of the minima of /. 
At low temperatures, the three maxima of the same 
height are obtained only for h = [Fig. GJa)]. They 
correspond to a triply degenerate ordered phase of the 
system for h = 0. By increasing temperature the max- 
imum at the origin (Fig. [3] b), which corresponds to 
the disordered state, appears and starts to rise. For 



Ao(3) = 4 In 2 = 2.772..., all the four maxima reach 
the same height (Fig. @|, indicating the coexistence be- 
tween disordered and ordered phases at the first-order 
transition point (K, h) = (Ao(3), 0). At still higher tem- 
peratures, K < Ao(3), the maxima behave differently 
for positive and negative values of h (Fig. [5|). For pos- 
itive values of the field, two maxima appear [Fig. EKb)], 
which represent the transition that we denote here as the 
transition, while for negative values of the field, 
there are three maxima [Fig. [5ja)[ representing the tran- 
sition denoted here as Isingj transition. 

At the KMS transition, the distance between the two 
maxima diminishes as K goes to its tricritical value, until 
they eventually merge at (Kg MS , hf r MS ) given by (23|) 
[Fig. Mb)}. 

The transition denoted as Ising, exists for all Ao(3) > 
K > K, 



Ising 



K (2) = 2. At A/ r smff and h 
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all three maxima [similar to those from the Fig. (B](a)], 
merge and the system undergoes a second-order phase 
transition of an Ising universality class. 

For a summary of the critical behavior of the present 
model let us turn back to the phase diagram in Fig. [5] 
There are four distinct phases: 

(I) Disordered phase at high temperatures, 

(II) For h > 0, there is an ordered phase in the Potts 
state 1. Ground state of this phase is non degenerate. 
Transition between this phase and disordered phase is of 
first order (the thick full line in Fig. [5]), except at the end 
point (tricritical point), where the transition is of second 
order, and where it is possible to go continuously from 
the phase I to the phase II. 

(III) For h — 0, there is a line of first-order transitions 
between the phase II and the phase IV (broken line in 
Fig. [2]). Its ground state is triply degenerate. 

(IV) For h < 0, the number of particles in state 1, cou- 
pled to the field, is suppressed favoring the particles in 
the Potts states 2 or 3. Ground state of this phase is 
twice degenerate. The first-order transition line separat- 
ing the phase IV and the disordered phase is marked by 
the open circles in Fig. O This line has no endpoint 
and goes down to h = — oo. The continuous transition 
between phases IV and the disordered phase is not pos- 
sible, except at infinity. 

Let us discuss first the scaling properties at different 
transitions. 



A. h > 0, KMS transitions 

A set of solutions satisfying the extrema conditions 
and (f3"01 corresponds to x = 0, while S is the solution 



1 After the approximation of Kihara et at, as explained in Sec. 
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FIG. 3. (Color online) The free-energy density -f(x,S) at (a) (K,h) = (3.5,0) and (b) (K, h) = (2.79,0). 




FIG. 4. (Color online) The free-energy density —f(x, S) at 
inverse temperature K — Kq(3) = 4 In 2 and h = 0. 



As the temperature rises, K < Kq(3), the distance be- 
tween the two maxima S± decreases and eventually van- 
ishes at the tricritical point 
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(endpoint of the thick full line in Fig. [5]). 

In the first-order transition regime and for large values 
of N, the dominant contribution to the partition function 
(|27p comes from the two minima of equal depth 
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Within this approximation, the calculation of zeros in the 
above expression, reduces to the following set of equa- 
tions for the free-energy densities at S± 



of the equation 
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The same equation is found in the approximation used by 
Kihara et al. [2ll. |27| to the second order phase transition 
in the MF Potts model, which neglects the fluctuations 
within the remaining [q — 1) states. For h and K, related 
through the Eq. J22]) (the thick full line in Fig. [2]) 
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the positions of the extrema S± [Fig. \5j[b)] ar e of the 
form 



S± = \ ± A S, 
where A S is a solution of the equation 
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(34) 



Re f{S-)-f(S+) =0, 



(37) 
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f{S-)-f(S+) =(2n+l)7r, n = 0,l,- 



To the leading order in 1/N, the solutions of the above 
equations are of the form 



1 4 N 8 V ' 



(38) 
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Remark that, to the leading order in 1/N, the real part 
of zeros does not depend on n. Numerical calculations 
(Figs. [7|and[8]) suggest that this remains true to all orders 
in 1/iV, i.e. that, at KMS transition, all the zeros lie on 
a straight line parallel to the imaginary axis. 

Calculation of the density of zeros at the Yang-Lee edge 
is performed in the same way as for the Ising (q = 2) case. 
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FIG. 5. (Color online) The free-energy density —f(x, S) at K 
-0.04169621964961 • • • and (b) Hkms = 0.0181471805599452 
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FIG. 6. (Color online) The free-energy density —f(x, S) at inverse temperature K 
(a) h = -0.06333826060 • • • and (b) h = h? r MS = 0.02648051389 • • • . 
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8/3 = 2.6 and external field: 



Relation ([16]) defines the density of zeros, while its value 
at the edge is given by the distance between the two zeros 
closest to the real axis. By inserting expansions f|38[) into 
(fTT)) . one gets 



2AS snf 1 
fjN = —. + O ( — 



3 71 



N 



(39) 



On the entire first-order transition line, A S > making 
g finite. By approaching the tricritical point, A S — > 0, 
and vanishing of g indicates the second-order transition 
at the tricritical point. 



B. h < 0, Ising transitions 

Negative values of the field bring the particles in the 
Potts state 1 (as defined in Eq. (|19|) ) to an energetically 
higher level than the particles in the states 2 and 3. Con- 
sequently, the system prefers to have most of the particles 
in states 2 and 3. In the limit h — > — oo. transitions of 



particles into state 1 are completely forbidden and the 
model reduces to the pure two-state (Ising) model with- 
out an external field. Thus, the negative field acts as a 
chemical potential: it regulates the number of particles 
in states 2 and 3. 

In the limit h — >• — oo the extrema conditions ([2TJ1 13X1)) 
and the free-energy density (|2"5)h reduce to 



•K 



1 



•K 



(40) 



-/(x,-i/2) = -jq- 



-ln 



2 cosh 



xK 



-const. (41) 



Equation P0|) has x — > solution at K = K^ m9 = 
Kq(2) = 2, so that the position of Ising tricritical point 
is 



raising 1 Ising 



(2, -co) 



(42) 



In the range of inverse temperatures Kq(2) < K < Kq(3) 
(open circles in Fig. [5]), the free-energy density has a 
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shape similar to the one shown in Fig. (Ha) , and the tran- 
sition is of first order. By approaching K K^ mg = 2, 
the first order character of the transition becomes weaker 
and, at K = 2, all three extrema of the free energy merge, 
and the transition changes its character into a second- 
order one. 

By studying the behavior of the free-energy density 
around the two edges of the Ising line of transitions, it is 
possible to calculate the shape of h = h(K) in the two 
limits. Close to the point (K, h) = (Kq(3), 0) 

h(K) ~ In ( e % - 1 V (43) 

Close to the point (K, h) = (Kq(2),— oo). 

h(K)~ hi (#~ 2 ) + y ( 44 ) 

In reference [27]], authors discuss the cases with h < 
and q > 4 proving the existence of the transition in the 
range 

K G (K (q-l),K (q)J , he (- oo,o). (45) 

Our results show explicitly how the above statement is 
extended to the case q = 3. 

1. Numerical results 

The numerical calculation of Yang-Lee zeros starts 
from Eq. ([23]) 

N 

Z N = a n (K) z n = Re Z N +i Im Z N = 0. (46) 

Numerical solutions of Re Zn = [blue (dark gray) lines] 
and Im Zn = [orange (light gray) lines] with their 
intersections giving the positions of Yang-Lee zeros are 
illustrated in Figs. [7] and [5] for two typical values of K. 
Two lines of zeros, similar to those in Fig. were also 
observed in a numerical study @ , of a three-dimensional 
q = 25 Potts model with nearest-neighbour interactions 
on small lattices. 

Following the same procedure as the one presented in 
Sec. [TTl we analyzed the positions of the two zeros clos- 
est to the real axis for a sequence of large system sizes 
ranging from N = 10 3 up to 10 6 sites. The position of 
the edge of the line of zeros and the edge density of ze- 
ros were extrapolated to the limit N — > oo by using the 
simple least-squares fit to the form 

const const 

hl,2,N = "l,2,extr + "T^J , §N = ffoxtr + -J-j— ■ (47) 

The results of these extrapolations in the three different 
regimes, corresponding to (K < K n (3),h > 0),(K > 
Kq(3), h = 0), and (K < Kq(3), h < 0), are presented in 



Tables HVl El and ED The K < K (3) regime contains 
two lines of transitions and in addition the Yang-Lee edge 
singularity for sufficiently high temperatures so that the 
gap around the real axis is oper0. At lower temperatures, 
K > Kq{2>), there is only one, h = 0, line of zeros. 

Table IIVI contains data describing the h > or KMS 
line of transitions (thick full line in Fig. [5]). 
At K > K^ MS , the imaginary part of the field, is equal 
to zero, which means that the gap is closed at low temper- 
atures. At higher temperatures, K < K^ MS , the zeros 
accumulate around the point (hi , hi > 0) which means 
that the gap is open with the Yang-Lee edge singularity 
at its edge. Real part of the edge still satisfies Eq. ([32| . 
Density of zeros at K > K^, MS has a finite value, and 
the transition is of first order. This finite value can be 
compared to the analytical expression (|39[) with high ac- 
curacy. At higher temperatures, K < K^ MS , the edge 
density of zeros vanishes and the transition is of second 
order. 

The set of convergence exponents y2,extr and y g ^xtr has 
the same values as the corresponding exponents of the 
q = 2 model (Table fl} , showing that it belongs to the 
same universality class. Since within the approximation 
by Kihara et a/.the handling of fluctuations in the model 
with q > 2 is essentially reduced to the calculation of 
fluctuations in the two-state model, this result is not sur- 
prising. 

One may clearly distinguish the tricritical exponents for 
K = K tr from the Yang-Lee edge exponents obtained for 
higher temperatures (K < Ktr)- The convergence expo- 
nent yi^extr is very close to 1, which is anticipated by Eq. 
([38| for the first order transitions. At K < K% MS , the 
transition is of second order, and we have no analytical 
expression for yi, but it seems that the value 1 remains 
for all temperatures. 

At low temperatures, K > -Ko(3) (see Table W\ . zeros 
of the partition function are the result of the competition 
between the three (at K > Ko(3), Fig. [3]) or the four (at 
K = A'o(3), Fig. Q| extrema of /. The transition is of 
a strong first-order type and loci of zeros form a straight 
line hi = (Figs. [5] a and[H]b). The convergence ex- 
ponents are equal to 1, except for the exponent of the 
edge density of zeros at Ao(3). A qualitatively different 
transition, discussed in Sec. MI B I appears for h < (Ta- 
ble IVI1 open circles in Fig. [5]) . The numerical precision 
of the results in Table I VII is much poorer than that in 
Table IIVI The main reason for that could be attributed 
to the fact that the corresponding zeros are fewer and 
more spaced so that even the closest two zeros are much 
more distant from the real axis, than they are for another 
line. For this reason, we do not proceed here with the 
calculation of the density of zeros used for other cases. 
The numerical extrapolations /ii,2,extr and /ii,2, exact pre- 
sented in Table I VII yield the exact values [given by ([2!?]) 



In present work, we analyzed it only around h > tricritical 
point. 
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TABLE IV. The results of extrapolation for h > 0: real and imaginary parts of the closest zeros, hi,a, and the edge density 
of zeros g, followed by the corresponding convergence exponents. The error bars are estimated to be of order of the last cited 
digit. The numerical part of the errors is at least one order of magnitude smaller than the overall errors (including also the 
finite-size corrections). 


if 




h\ f extr 


^1, exact 


J/l,extr 


%fl , guess 


if o(3) = 2.772 ■ 




10" 12 





1.000003 


1 


2.7 




0.018147180556 


0.018147180559 • • • 


1.000003 


1 


iftr = 2.6 




0.026480513887 


0.026480513893 • • ■ 


1.000004 


1 


2.6 




0.04314719 


0.04314718 ■ ■ ■ 


0.99993 


1 


2.5 




0.0681470 


0.0681472 ■ • • 


1.0005 


1 


2.4 




0.0931470 


0.0931472 • • • 


1.0006 


1 


2.3 




0.1181471 


0.1181472 ■ ■ ■ 


1.0005 


1 


if 




ft-2,extr 


h<2, exact 


J/2,oxtr 


J/2, guess 


if o(3) = 2.772 ■ 




10" 8 





1.001 


1 


2.7 




10 





1.01 


1 


T JC ft f Q „ ' 

Ktr = 2.6 




10" 8 





0.749 


3/4 


2.6 




0.00533 




0.668 


2/3 


2.5 




0.02126 




0.663 


2/3 


2.4 




0.04347 




0.658 


2/3 


2.3 




0.07091 




0.655 


2/3 


if 




ffextr 


Pexact 


J/g,cxtr 


Va, guess 


if o(3) = 2.772- 




0.05305168 


0.05305164 • • • 


0.997 


1 


2.7 




0.03048 


0.03047 ■ • • 


0.94 


1 


Kf r MS = 2.6 




HT 6 





0.2503 


1/4 


2.6 




10 -3 





0.28 


1/3 


2.5 




10" 3 





0.29 


1/3 


2.4 




10 -3 





0.335 


1/3 


2.3 




10~ 3 





0.29 


1/3 


TABLE V. The results of extrapolation for if > ifo(3), h = 0: the real and imaginary part of the closest zeros, /ii,2 and the 
edge density of zeros g, followed by the corresponding convergence exponents. The error bars are estimated to be of order of 
the last cited digit. The numerical part of errors is by at least one order of magnitude smaller than the overall errors (including 
also the finite-size corrections). 


if 




ftl.extr 


h\ ; exact 


?/l,extr 


yi .guess 


3.5 




10 " 7 





1.001 


1 


3.0 




10- 7 





1.007 


1 


ifo(3) = 4 In 2 = 


= 2.772 ■ • ■ 


10 ~ 5 





1.05 


1 


if 




/l2,extr 


^2, exact 


?/2,extr 


]J2, guess 


3.5 




10 " 7 





1.0007 


1 


3.0 




10 -° 





1.004 


1 


ifo(3) = 4 In 2 = 


= 2.772 ■ • ■ 


10 ~ 5 





1.04 


1 


if 




Qext r 


Pexact 


,extr 


Vg. guess 


3.5 




0.13836168 




0.997 


1 


3.0 




0.11401535 




0.992 


1 


if o(3) = 4 In 2 = 


= 2.772 ■ ■ ■ 


0.08169236 




1.53 


3/2 



10 




-0.1 -0.05 0.05 -0.05 -0.04 -0.03 -0.02 -0.01 0.01 0.02 



h l h l 

FIG. 7. (Color online) q — 3 and N = 3 000. The positions of the part of the Yang- Lee zeros closest to the real field axis are 
given as the intersections of the blue (dark gray) and orange (light gray) lines: (a) K — 2.6 and (b) K — 2.7. The broken black 
line is the KMS value of the real field predicted by relation (|32l) . 




FIG. 8. (Color online) q — 3 and N = 10 000. The positions of the part of Yang- Lee zeros closest to the real field axis are given 
as the intersections of the blue (dark gray) and orange (light gray) lines: (a) K — 7^0(3) = 2.772 . . . and (b) K — 3.5. 



TABLE VI. The results of extrapolation for h < 0: the real and imaginary part of the closest zeros, hi t 2, followed by the 
corresponding convergence exponents. The error bars are estimated to be of order of the last cited digit. The numerical part 
of errors is by at least one order of magnitude smaller than the overall errors (including also the finite-size corrections). 



K 


^l,extr 


/^l, exact 


yi.cxtr 


2.7 


-0.04162 


-0.04169 ■ ■ ■ 


1.205 


2.6 


-0.1093 


-0.1099' •• 


0.685 


2.5 


-0.192 


-0.193- ■■ 


0.479 


K 


^2,extr 


^2, exact 


2/2,cxtr 


2.7 


10" 4 





1.068 


2.6 


HT 4 





0.79 


2.5 


10 -3 





0.69 
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and (|3T)|) ] with numerical precision up to three digits 
or more. The convergence exponents yi,2,extr differ, 
however, significantly from the scaling exponent of the 
first-order phase transition, which should persist up to 
h = — oo, as argued in Sec. MI Bl 



IV. CONCLUSION 

The phase diagram of the two- and three-state Potts 
model with infinite-range interactions in the external field 
was analyzed by studying the partition function zeros in 
the complex field plane. 

In the three-state case, we derived the exact two-pa- 
rameter expression for the free-energy density after a 
twofold application of the Hubbard-Stratonovich trans- 
formation. By applying the saddle-point approxima- 
tion in the two-parameter space, we reproduced well the 
known analytical results in different regimes of the phase 
diagram with zero and positive external field. Calcula- 
tions could also be extended to the negative fields, where 
another, Ising-like, phase transition occurs. We derived 
an implicit analytical expression and performed exact nu- 
merical calculation for the critical line of this Ising-like 
transition, showing that it is of first order, and becomes 
of second order only in the limit h — > — oo, unlike the case 
of the three-dimensional Potts model with short-range in- 
teractions, where the tricritical value of the field is finite 




The same procedure was efficient for the study of the 
partition function zeros in the complex field plane. We 
have shown, that in this model, with a rather complex 
phase diagram and with consequently more complicated 



loci of zeros including multiple lines, it is still possible 
to identify different regimes of critical behavior from the 
size scaling properties of the few zeros closest to the real 
axis. We successfully identified the regime of the first- 
order phase transitions, and well reproduced the critical 
exponents belonging to the two second-order phase tran- 
sitions, the tricritical point, and the Yang-Lee edge sin- 
gularity. In the case of the tricritical point, an unusual 
feature is that the line of zeros reaches the real axis at 
the point with nonzero field (i.e. loci of zeros lie on the 
non unit circle in the complex activity plane) . The Yang- 
Lee edge singularity could be observed for temperatures 
above the tricritical point temperature (K < K tr ), and 
was found to belong to the same universality class as that 
of the MF Ising model. 

Due to the fact that the external field h breaks the 
symmetry between Potts states of the q = 3 model in 
essentially the same way as for all q > 3 models, one can 
expect that the general shape of the phase diagram pre- 
sented in figure [5] as well as a major part of the obtained 
results could be extended to models with q > 3. 

It would also be interesting to perform a complemen- 
tary study of a more general case of the Yang-Lee zeros 
for the Potts model with infinite range interactions with 
the power-law decay, where some previous work has been 
performed [3(| for the Ising model by the finite-range 
scaling approach (3l| . 
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